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Using the time-dependent mean-field Gross-Pitaevskii equation we study the dynamics of small 
repulsive Bose-Einstein condensed vortex states of 85 Rb atoms in a cylindrical trap with low angular 
momentum fiL per atom (L < 6), when the atomic interaction is suddenly turned attractive by 
manipulating the external magnetic field near a Feshbach resonance. Consequently, the condensate 
collapses and ejects atoms via explosion and a remnant condensate with a smaller number of atoms 
emerges that survives for a long time. Detail of this collapse and explosion is compared critically 
with a similar experiment performed with zero angular momentum (L = 0). Suggestion for future 
experiment with vortex state is made. 



PACS numbers: 03.75.Fi 
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I. INTRODUCTION 

Recent observation [|, ^ of Bose-Einstein condensates 
(BEC) of dilute trapped bosonic atoms with repulsive 
interaction at ultra-low temperature has intensified the- 
oretical and experimental studies on various aspects of 
the condensate S. One fascinating feature is the obser- 
vation of quantized vortices || and vortex lattice (5[ ||] 
in the condensate as this is intrinsically related to the 
existence of superfluidity. Another interesting feature is 
the formation of a stable condensate composed of a finite 
number of attractive atoms less than a critical number 
Acr • The third noteworthy feature is the observation 
of Feshbach resonances in 23 Na §, 85 Rb § and Cs § 
atoms, as in the presence of such resonances the effective 
atomic interaction can be varied in a controlled fashion 
by an external (background) magnetic field jio) . 

For superfluid 4 He(II) in a rotating container, no mo- 
tion of the fluid is observed below a critical rotational 
frequency. Above this frequency quantized vortices ap- 
pear in 4 He(II) manifesting its superfluidity. However, 
because of the strong interaction, the theoretical de- 
scription of this system is not easy. Quantized vor- 
tices have been observed ||, ||, [| in trapped BEC 
and can be generated in theoretical mean-field models 
P, PLP, 0= 0, based on the Gross-Pitaevskii 
(GP) |1 gflequation. Different ways for generating vortices 
in a BEC have been suggested [O, e.g., via spontaneous 
formation in evaporative cooling |14| , via controlled exci- 
tation to an excited state JT|| , by stirring a BEC using a 
laser with an angular frequency above a critical value p2[ , 
or by the rotation of an axially symmetric trap with an 
angular frequency above a similar critical value |[(|. In 
contrast to liquid 4 He(H), a trapped BEC of small size is 
dilute and weakly interacting, which makes a mean-field 
analysis appropriate. 

The observation of a condensate of attractive 7 Li atoms 
and the subsequent measurement of the critical number 
Acr is in good agreement with the mean field analyses 
in a spherically symmetric trap J19| , although the agree- 



ment is not as good in the case of 85 Rb [^o) atoms in an 
axially symmetric trap . If the number of atoms can 
somehow be increased beyond this critical number, due 
to interatomic attraction the condensate collapses emit- 
ting atoms until the number of atoms is reduced below 
Acr- With a supply of atoms from an external source 
the condensate grows again beyond Acr and a sequence 
of collapses has been observed in 7 Li by Gerton et al at 
Rice |g , where the number of atoms remain close to Acr 
and the collapse is driven by a stochastic process. 

Recently, a more challenging experiment was per- 
formed by Donley et al. @ at JILA on a BEC of 85 Rb 
atoms [|| in an axially symmetric trap, where they var- 
ied the interatomic interaction by an external magnetic 
field near a Feshbach resonance (To| . Consequently, they 
were able to change the sign of the atomic scattering 
length, thus transforming a repulsive condensate of 85 Rb 
atoms into a collapsing and highly explosive attractive 
condensate and studied the dynamics of the same p2[ . 
Immediately after the jump in the scattering length, one 
has a highly unstable BEC, where the number of atoms 
could be much larger than Acr- Donley et al. |Q have 
provided a quantitative estimate of the explosion by mea- 
suring the number of atoms remaining in the BEC as a 
function of time until an equilibrium is reached. Because 
this phenomenon of emission of a very large number of 
atoms in a small interval of time is reminiscent of an ex- 
plosion and looks very much like a tiny supernova, or 
exploding star, the researchers dubbed it a "Bosenova" . 

The essential aspects of the above experiments at Rice 
by Gerton et al [|| and at JILA by Donley et al have 
been theoretically described by a variety of authors us- 
ing the GP equation [jl8|. The theoretical analyses have 
not only produced the time-independent results, such as, 
the critical number Acr 119, 21 , but also time-dependent 



results, such as, the variation of number of atoms of the 
BEC during collapse and explosion Q ||[ ||, ||, |27], ||) , 
both in reasonable agreement with experiment. This con- 
solidates the use of the mean-field GP equation in de- 
scribing the dynamics of collapsing and exploding BEC 
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of small to medium size. These BEC's composed of sev- 
eral thousand atoms can be considered dilute and weakly 
interacting and hence amenable to mean-field treatment. 
Motivated by the above success, using the GP equation 
we propose the numerical simulation of the dynamics of 
a rotating Bosenova with a single vortex composed of a 
small number (several thousands) of 85 Rb atoms as in 
the experiment at JILA |g2| . We consider a single vortex 
state |ll], |2l| , as appropriate for small condensates, op- 
posed to vortex lattice for large condensates || [| |l7J . A 
comparison of the present results with future experiment 
will provide a more stringent test for the mean-field GP 
equation. 



In this paper we perform a mean-field analysis based 
on the time-dependent GP equation to understand the 
collapse and explosion of the attractive vortex state of 
85 Rb atoms in an axially symmetric trap. To account 
for the loss of atoms from the strongly attractive con- 
densate we include an absorptive nonlinear three-body 
recombination term in the GP equation. The three-body 
recombination rate we use in numerical simulation is the 
same as used in a similar study with BEC's with zero an- 
gular momentum [ p7| and is in agreement with previous 
experimental measurement [^9| and theoretical calcula- 
In the present investigation we consider the 
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complete numerical solution of the mean-field GP equa- 
tion for an axially symmetric trap as in the experiment 
at JILA ]2^| . As in the experiment at JILA with non- 
rotating condensates, we find that, also in the case of 
rotating vortex states, a large number of atoms could be 
emitted in a small interval of time and one could have a 
Bosenova-type explosion. 



Throughout the present numerical simulation we make 
the assumption that the axial symmetry of the system is 
maintained. For small values of nonlinearity a dynamical 
quadrupole instability may cause an attractive BEC vor- 
tex state to split into two pieces that rotate around the 
axial direction 31 . These pieces may unite to recover the 
original vortex and this split-merge cycle repeats. Simi- 



lar instability is known to exist for BEC in a toroidal trap 
p2| . Clearly, a full three-dimensional calculation of the 
collapsing phenomena taking into consideration the effect 
of the splitting of vortex states seems to be practically 
impossible at present and will be a welcome future work. 
In view of this here we present an axisymmetric model 
of the same which is expected to provide the essentials of 
the collapse dynamics of the vortex states. 



In Sec. II we present the theoretical model and the nu- 
merical method for its solution. In Sec, III we present our 
results. Finally, in Sec. IV we present a brief discussion 
and concluding remarks. 



II. THE GROSS-PITAEVSKII EQUATION 

A. Theoretical Model 

The time-dependent Bose-Einstein condensate wave 
function \&(r;r) at position r and time r allowing for 
atomic loss may be described by the following mean-field 
nonlinear GP equation 0, lill 
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K 3 N 2 \^(r;r)\ 4 ) 



*(r;r)=0. (2.1) 



Here m is the mass and N the number of atoms in the 
condensate, Q = 4Trh 2 a/m the strength of interatomic in- 
teraction, with a the atomic scattering length. The terms 
K2 and K3 denote two-body dipolar and three-body re- 
combination loss-rate coefficients, respectively. There are 
many ways to account for the loss mechanism (23[ ^4| . 
Here we simulate the atom loss via the most important 
quintic three-body term K3 |2^, ^6|. The contribu- 
tion of the cubic two-body loss term ]29| is expected to 
be negligible |2^, ||(| compared to the three-body term 
in the present problem of the collapsed condensate with 
large density and will not be considered here. 

The trap potential with cylindrical symmetry may be 
written as V(r) = ^mui 2 (r 2 + \ 2 z 2 ) where u> is the an- 
gular frequency in the radial direction r and \u> that in 
the axial direction z, with A the aspect ratio. We are 
using the cylindrical coordinate system r = (r, 9, z) with 
9 the azimuthal angle. The normalization condition of 
the wave function is J dr|^(r;r)| 2 = 1. 

The GP equation can easily accommodate quantized 
vortex states with rotational motion of the BEC around 
the z axis. In such a vortex the atoms flow with tan- 
gential velocity Lh/(mr) such that each atom has quan- 
tized angular momentum Lh along the z axis. This cor- 
responds to an angular dependence of 



*(r, t) = ip(r, z, t) exp(iL9) 



(2.2) 



of the wave function, where exp(iL9) are the circular 
harmonics in two dimensions. 

Now transforming to dimensionless v ariables defined 
by x = \/2r/l, y = y2z/l, t = tuj, I = yjh/{muj), and 
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we get from Eqs. (gTj) and (gj) 
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<p{x,y,t)=0, (2.4) 



where n = Na/l, k = 8V2%, and £ = 4K 3 /(a 2 l 4 uj). From 
theoretical p3| and experimental p9| studies it has been 
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found that for negative a, K 3 increases rapidly as \a\ n , 
where the theoretical study favors n — 2 and we represent 
this variation via this quadratic dependence. This makes 
the parameter £ above a constant |B7[ for an experimental 
set up with fixed I and u> and in the present study we 
employ a constant £. 

The normalization condition of the wave function for 
K 3 = is 



■A/norm = 27r 



" "° dx f°° 

,no]-!i) -ri — / dy\tp(x,y;t)\ 2 = 1. (2.5) 

'0 ^ J-oo 



However, in the presence of loss K 3 > 0, A/norm < 1. 
The number of remaining atoms N in the condensate is 
given by N = A/norm > where Nq is the initial number. 

The root mean square (rms) sizes Xrms an d 2/rms are 
defined by 



2tt 



oo poo 



x rms — »/■ 

yvnorm Jo ./-oo 

2 _ 2tt /-""da; ^ 
2/rms — 



/ dy|<£>(a:,y;i)| a, (2.6) 



norm 



/ 2£ / dj,|y,(x jy ;t)|V. (2-7) 

JO x J-oo 



B. Calculational Detail 



We solve the GP equation (2.4) numerically using 
a split-step timc-itcration method using the Crank- 



Nicholson discretization scheme J21], |27j, |34|, |35J, |36| . We 

discretize the GP equation with time step 0.001 and space 
step 0.1 spanning x from to 15 and y from —35 to 35. 

It is now appropriate to calculate t he pa rameters of the 
present dimensionless GP equation (2.4) corresponding 
to the experiment at JILA for L = [22[. As in that 
experiment we take the radial and axial trap frequencies 
to be ^ rac ji a i = 17.5 Hz and ^ ax j a j = 6.8 Hz, respectively, 
leading to A = 0.389. The harmonic oscillator length I of 
85 Rb ato ms for u j = 2?r x 17.5 Hz and m m 79176 MeV 
is I = y/h/(muj) = 26070 A. One unit of time t of Eq. 
Q is 1/w or 0.009095 s. 

We consider a stable 85 Rb condensate of No = 16000 
atoms with scattering length a m jti a l = 7a , a = 0.5292 
A. This wave function is obtained by time iteration of Eq. 
( |2.4D employing the following normalized initial solution 
with a single central vortex pjj 



<p(x,y) 



A 



-,1/4 



2 2L + 3 tt 3 (|L|!) 2 



x l+\L\ e -(x- , +Xy")/4 



(2.8) 



for n = 0. In the course of above time iteration the 
nonlinear parameter n was increased by steps of 0.0001 
until its final value is attained. Then during an interval 
of time 0.1 ms the scattering length was ramped to a = 
a collapse" absorptive term £ was set equal to zero 
throughout above time iteration. 

The final condensate is strongly attractive and un- 
stable and undergoes a sequence of collapse and ex- 
plosion. In our numerical simulation with L =/= 



we consider a set of different values of ^collapse ( = 
— 263a , — 100ao, — 30a , — 20ao, etc.) as well as N = 
6000. 

For the simulation of collapse and explosion a nonzero 
value of £ (= 2) is chosen for different a m iti a i, a collapse' 
and No as in Ref. and the time-evolution of the 

GP equation is continued. This value of £ reproduced 
the essentials of the experiment at JILA [B2| reasonably 
well for L = as well as produced J2lj a K3 in reasonable 
agreement with a previous experiment (K3 = 4.24x 10~ 25 
cm 6 /s) 29| and theoretical calculation (K 3 = 6.7 x 10~ 25 
cm 6 /s) [B0] for a = — 370a . In particular we use K 3 = 
9 x 10~ cm 6 /s for a — — 370ao. For smaller values 
of \a\, the K3 values are scaled down using the relation 
K 3 oc a 2 . 



III. NUMERICAL RESULT 



The numerical simulation using Eq. (2.4) with a 
nonzero £ as described above immediately yields the re- 
maining number of atoms in the condensate after the 
jump in scattering length. The remaining number of 
atoms vs. time is plotted in Fig. 1 (a) for o m j|.j a j = 7ao, 
a collapse = — 30ao and — 263ao, A*o = 16000, and 
L = 0, 1,2,4, and 6. In Fig. 1 (b) the same results for 
A^o = 6000 are plotted. In this figure we also plot some 
results of experiment at JILA for L = (2^, ^7|. These 
experimental results are in agreement with the simulation 
for L = 0. In Fig. 2 (a) we plot the particle loss curves 
for A^o = 6000, aj n ^i a l = 7ao an d different & C0 H a p Se for 
L = 0. In Fig. 2 (b) we plot the same for L = 1. 

From Figs. 1 and 2 we find that for a fixed TVo, for 
a sufficiently small L or a sufficiently large I ^collapse I ' 
there could be collapse and explosion during a relatively 
short interval of time (called decay time) with the loss of a 
large fraction of the atoms. However, there is no collapse 
for a large enough L or a small enough I a collapse I ' ^ or 
example, for A*o = 6000 in Fig. 1 (b) there is no collapse 
for L > 2 for |a co n a p Se | = 30ao and in Fig. 2 (a) there 
is no collapse for |» co n a p Se | < 5ao for L = 0. 

In the experiment at JILA |2^] for L = it was 
observed that the strongly attractive condensate after 
preparation remains stable with a constant number of 
atoms for an interval of time i C0 Hapse' caue d collapse 
time. This behavior is physically expected for medium 
to small values of I ^collapse I 50°o)- Immediately af- 
ter the jump in the scattering length to a negative value, 
the attractive condensate shrinks in size during £ C0 Hapse' 
until the central density increases to a maximum. Then 
the absorptive three-body term takes full control to ini- 
tiate the explosion that last for few milliseconds. Conse- 
quently, the number of atoms remains constant for time 
* < ^collapse' present results in Figs. 1(a) and (b) 
also show this behavior for I a collapse I = 30ao- However, 
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FIG. 1: The number of remaining atoms N(t) in the conden- 
sate of (a) No = 16000 and (b) 6000 atoms after ramping the 
scattering length from a^tia] = 7a to a co llapse ~ ™ 30ao 
(dashed line) and — 263ao (full line) for different I as a func- 
tion of time t. The curves are labeled by their respective L 
values. Solid circles represent results of experiment at JILA 
Q for L — 0, No = 16000 and a coYl e = -30a and 
dashed-dotted lines represent the average [E7J over experimen- 
tal results at JILA j22j for L = and Q co n a p Se = — 263ao 

for larger I ^collapse I ( = 263ao), the atomic attraction is 
very strong and the central density increases to a maxi- 
mum quickly to start the explosion and * co llapse * s c l° se 
to zero. In Figs. 2 (a) and (b) we see the dependence of 
particle loss and i co n a p Se on |« C0 H a p Se l f° r = 6000, 
a initial = ^ ao > anc ^ L = and L = 1, respectively. From 



FIG. 2: The number of remaining atoms N(t) in the conden- 
sate of No = 6000 atoms for (a) L = and (b) L = 1 after 
ramping the scattering length from Gtj n jtial = ^° different 
final o co n a p Se as a function of time t. The curves are labeled 
by their respective a co n a p Se values. 

Figs. 1 and 2 we find that i co llapse mcreases with L for 
a fixed a co n apS e and with a collapse for a fixed L. 

From Figs. 1 and 2 we find that after the collapse, the 
number of particles drop sharply during a small interval 
of time called decay time (few milliseconds) , which means 
that the condensate emits a large number particles in an 
explosive fashion. This emission of particles is termed 
explosion. 

After a sequence of collapse and explosion, for L = 
Donley et al. |2^] observed a "remnant" condensate 
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FIG. 3: The number of remaining atoms N(t) in the con- 
densate of iVo = 16000 atoms for different L after ramp- 
ing the scattering length from aj n jtial ~ ^° ^ ne fi na -l 
"collapse = ~30ao as a function of time t. The curves are 
labeled by their respective L values. 



FIG. 4: The collapse time i co llapse vs ' l a collapse l/ a ° f° r 
"initial = ^ anc ^ ~ 6000 for different L. Solid circle with 
error bar: experiment Jz| for L — 0; open circle: axially 
symmetric mean-field model of Ref. [E6| for L — 0; arrows 
are the acr/ao values, full line: present theory for different L. 



of -^remnant atoms at large times containing a fraction 
of the initial No atoms. Figures 1 and 2 show such a 
behavior for different values of L and a C ollapse' In all 
cases the decay time during which the explosion takes 
place is small and of the order of few milliseconds. The 
decay time for vortex states (L ^ 0) is smaller than for 
nonrotating condensates (L = 0). 

We studied the time evolution of the condensate for 
larger times. In Fig. 3 we plot the loss curves for 
A = 16000, a initial = 7a , a C ollapse = - 30a o> and 
L = 0,1,2,4,6 at larger times. The BEC continues to 
lose atoms if left for a long time but at a rate much 
slower than during the first explosion that we call pri- 
mary. However, in the experiment at JILA Donley et al 
observed that a remnant condensate containing a fraction 
of the atoms survived with nearly constant number for 
more than one second. One possible reason for this dis- 
crepancy could be the following. In the actual experiment 
at JILA a major portion of the emitted atoms, called the 
burst atoms, remain trapped and oscillate around the 
central remnant. The presence of the burst atoms make 
the measurement of the number of atoms in the rem- 
nant a difficult task |37j . Some of these burst atoms may 
also rejoin the remnant to compensate for the three-body 
loss at large times. Such an effect is not included in the 
present model, which, hence, presents a larger loss for 
the remnant compared to experiment. 

We also observe an interesting phenomenon in Fig. 3, 
e. g., the occurrence of smaller secondary and tertiary ex- 
plosions after the primary one observed for small times. 



For L = 2 after the primary collapse and explosion at 
t < 10 ms with the loss of about 7000 atoms, there is 
another collapse and explosion with loss of about 2000 
atoms at f w 150 ms. The primary and secondary ex- 
plosions are separated by a large interval of time. We 
also see much weaker explosion(s) in the course of time 
evolution in Fig. 3, where the particle number varies in 
small steps. These explosions could be termed tertiary 
with the loss of few hundred atoms. It might be inter- 
esting to see if such secondary and tertiary explosion(s) 
could be observed experimentally. 

Donley et al. |22| provided a quantitative measure- 
ment of the variation of collapse time i co llapse w ^h the 
final scattering length « C ollapse for a given fflj^ial = 0, 
Ao = 6000, and L = 0. We calculated this variation 
using our model for L = 0, 1 and 2. The ^collapse vs ' 
l a collapsel/ a ° pl°t s f° r L = 0,1,2 are exhibited in Fig. 
4 and compared with experimental data for L = |2^] 
as well as with another calculation using the mean-field 
GP equation in an axially symmetric trap for L — p6[ . 
We see i co n a p Se decreases with l a co ll a p Se l/ a o starting 
from an infinite value at |a co n a p Se | = acr, where acr 
is the minimum value of l a co n a p Se l that leads to col- 
lapse and explosion. The critical value acr increases with 
L and so does ^collapse ^ or a fi xec l l a collapsel' -^ or a 
given Ao, a critical value of n = ncr can be defined via 
ncr = Noaci/l- A necessary condition for collapse is 
Ao|a co rj a p Se |/Z > ricr |H|- The value of ncr for a specific 
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case is calculated as in Ref. [Elj. Consequently, acr/ao 
can be obtained. The value of acr/ a o so evaluated for a 
specific L is shown by an arrow near the curve for that 
particular L in Fig. 4. The i collapse vs. |a co ll a psel/ a o 
curves should tend to infinity at the respective arrows 
and they do so in Fig. 4. There should not be any col- 
lapse for |a collapse | < a CT . 
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FIG. 5: Remnant number vs. initial number for a j n jtia,l = 
7ao, L — 0, 1 and 2, and different a co n a p Se - The results 
are represented by different types of triangle, circle, square, 
and inverted triangle for L = 0, 1, and 2 and a co ]l a p Se — 
— 21ao, — 30ao, — lOOoo and — 255ao as indicated in the figure. 



Donley et al. (2 
atoms for L = 



^collapse 



measured the number of remnant 
a initial = an< ^ different Nq and 
We plot the same in Fig. 5 for L — 0, 1 

and 2. The remnant number plotted in this figure is the 
number after the primary explosion(s) and not during or 
after possible secondary and tertiary explosions at larger 
times. In Figs. I and 2 the remnant number is obtained 
around t ~ 20 — 30 ms and not at few hundred millisec- 
onds. Our results in Fig. 5 for L = agree well with 
the measurements of Donley et al. |22[ ]. In general the 
remnant number decreases with increasing L. However, 
there are some cases where the opposite trend has been 
observed in Figs. 1 and 2 and as well as in Fig. 5. For the 
smallest values of Nq in Fig. 5, the condensate remains 
stable for L > 0, and there is no collapse and explosion 
and hence no remnant numbers for L = I and 2. For 
certain iVo the results for all three L's are not plotted as 
they coincide with other remnant numbers. The remnant 
number in some cases could be much larger than A^cr for 
times on the order of tens of milliseconds. 

Donley et al. [^2| observed that for L = the rem- 
nant condensate always oscillated in a highly excited col- 
lective state with approximate frequencies 2i/ ax j a j and 
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FIG. 6: The rms sizes irms (full line) and j/rms (dashed line) 
after the jump in the scattering length of a BEC of 16000 85 Rb 
atoms for L — 2 from <3j n jtial = *° a collapse ~ — 30<xo as 
functions of time t. 



"^radial being predominantly excited. This behavior 
emerges from the present simulation for all values of L. 
To illustrate this we plot in Fig. 6 sizes irms and yrms 
vs. time for the condensate after the jump in the scatter- 
ing length to — 30ao from 7oo for Nq = 16000 and L = 2. 
We find a periodic oscillation in Xi-ms an d yrms with fre- 
quencies 13.6 Hz (~ 2f ax j a j) and 35 Hz (~ 2^ ra( jj a i), 
respectively, as observed in experiment. 



IV. CONCLUSION 

In conclusion, we have employed a numerical simula- 
tion based on the solution |2l[] of the mean-field Gross- 
Pitaevskii equation ]l8| with cylindrical symmetry to 
study the dynamics of collapse and explosion (2^, |tJ of 
small attractive vortex states with L > 0. The explosion 
is initiated by a sudden jump in the scattering length 
from a positive to negative value exploiting a Feshbach 
resonance 0, ||, [hJ. In the GP equation we include 
a quintic three-body nonlinear recombination loss term 
p3 , |2~i| , p5| , |26| that accounts for the decay of the strongly 
attractive condensate. The results of the present simula- 
tion are to be considered as an extension of the experi- 
ment at JILA for L — |22 to small vortex states. 



We find the following features of this dynamics from 
the present numerical simulation: (1) The condensate 
undergoes collapse and explosion during a small interval 
of time of few milliseconds and finally stabilizes to a rem- 
nant condensate containing a fraction of initial number 
of atoms. The number in the remnant condensate for 
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times on the order of tens of milliseconds can be much 
larger than the critical number for collapse N CT for the 
same atomic interaction. (2) In some cases after the pri- 
mary explosion small secondary and tertiary explosions 
are observed. (3) The explosion takes place during a 
decay time of few milliseconds. This decay time for ro- 
tating Bosenova with vortex (L > 0) is smaller than the 
same for nonrotating condensate with L = 0. (4) The 
remnant condensate executes radial and axial oscillations 
in a highly excited collective state for a long time with 

frequencies ^ v Y a,6id\ anc ^ ^axiak (^) After the sudden 
change in the scattering length to a large negative value, 
the condensate needs an interval of time i C0 Uapse before 
it experiences loss via explosion. The interval i C0 H a p Se 
increases with L and a co n a p Se - 

The simulation of the particle loss in strongly attrac- 



tive rotating Bosenova, with a single axial vortex of small 
angular momentum per particle, may stimulate further 
theoretical and experimental studies. We have consid- 
ered small vortex states as they can be well described by 
the mean-field GP equation. This will provide a test for 
the usefulness of this equation in handling particle loss. 
Otherwise, a similar study with a large Bose-Einstein 
condensed vortex lattice ^| is more challenging from 
both experimental and theoretical points. 
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